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We develop a new tool to generate statistically precise dark matter maps from the cosmic mag- 
nification of galaxies with distance estimates. We show how to overcome the intrinsic clustering 
problem using the slope of the luminosity function, because magnificability changes strongly over 
the luminosity function, while intrinsic clustering only changes weakly. This may allow precision 
cosmology beyond most current systematic limitations. SKA is able to reconstruct projected matter 
density map at smoothing scale ~ 10 with S/N> 1, at the rate of 200-4000 deg 2 per year, depending 
on the abundance and evolution of 21cm emitting galaxies. This power of mapping dark matter is 
comparable to, or even better than that of cosmic shear from deep optical surveys or 21cm surveys. 
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Introduction. — The precision mapping of the universe, 
and the accurate determination of cosmological parame- 
ters have been enabled by the recent generation of cos- 
mic microwave background(CMB) experiments, galaxy 
and lensing surveys, and new analysis techniques. Weak 
gravitational lensing has emerged with a promising fu- 
ture of mapping dark matter directly, which would allow 
the inference of the state of the universe, including its 
dynamics and the nature of dark energy. Lensing is free 
from modeling assumptions, and can be accurately pre- 
dicted from first principles. Several major surveys are 
underway, under construction or in the planning stage. 
Currently, most attention has focused on using the lens- 
ing induced cosmic shear 1]. But such an approach is 
subject to a series of difficult experimental systematics 
0. CMB lensing^^ and 21cm background lensing Q are 
promising. But contaminations such as the kinetic Sun- 
yaev Zeldovich effect 5] and/or non-Gaussianity may de- 
grade their accuracy. In this paper we will address an 
alternative approach, the lensing induced cosmic magni- 
fication, which is not subject to the known problems, and 
could provide a robust statistical signal. 

Traditionally, intrinsic clustering had presented a se- 
rious problem to measurement of cosmic magnification. 
The observable quantity is the surface density of galaxies 
above some flux threshold. A variation in this surface 
density is then interpreted as lensing. Unfortunately, in- 
trinsic clustering is usually larger than the lensing in- 
duced signal. By utilizing the redshift information, in- 
trinsic clustering can be effectively eliminated in lensing 
correlation functions0,0- In this paper, we further show 
that, beyond the above statistical lensing measurement, 
2D convergence k maps can be reconstructed with lower 
systematics and larger sky coverage than cosmic shear 
maps, by utilizing both the redshift and flux informa- 
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tion of galaxies. 2D k maps not only provide indepen- 
dent and robust constraints on cosmology, but also are 
complementary to traditional shear maps. It allows one 
to explicitly and locally solve for non-reduced shear, an 
independent mode of checking E-B decomposition, and 
break the mass-sheet degeneracy p. 

Cosmic magnification. — Cosmic magnification causes 
coherent changes in the apparent galaxy number density. 
Let Nij be the observed number of galaxies (including 
false peaks) at the i-th flux bin and j-th redshift bin, 
falling into an angular pixel centered at direction n with 
angular size 0. It can be expressed as 

Nijih) = N i:j + N r l3 [WijKjin) + S gtij (h)] + 5N P ^(h) .(1) 

The signal Wk has unique dependence on galaxy flux 
through W = 2(a- 1). Here, a = -dhi[dn/dF]/dlnF - 
1 and dn/dF is the mean number of observed galaxies 
per flux interval N tj = JV£ + N^, N[j, N f %j are the 
mean number of detections, real galaxies and false peaks, 
respectively. S g and 5Np are galaxy intrinsic clustering 
and Poisson fluctuation, respectively. 

Our goal is to recover k of each angular pixel, given 
observables jV f - ; -, Nf f . IT and A 7 ,-, 20] . We consider 
SKA [2l]], which can detect ~ 10 8 high z galaxies through 
the neutral hydrogen 21cm emission line, n has typical 
value ~ 1%. To beat down Poisson fluctuations, > 10 4 
galaxies per angular pixel are required. Traditionally, ob- 
jects are selected at a 5c cut, where one can neglect the 
fraction of false detections. This of course also discards 
the majority of the signal. With a 0.5-er cut, one can re- 
duce Poisson noise at ~ 10 . To increase lensing signal 
while reducing S g contamination, we focus on source red- 
shifts z > 2. After averaging over the full redshift range 
z > 2, 8 g is still several times larger than k. However, 
SNp and 8 g ^j have different flux dependence to that of 
the signal. Weighting each galaxies by some function of 
their flux can suppress the prefactors of S g and SNp. In- 
tuitively, Eq. 1 implies the optimal estimator to be linear 
in Nu. 
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The predictions rely on the assumed HI mass func- 
tion n(Mni,z). We extrapolate the locally observed 

n(M m ,z) = n (z)(M m /M i ,)- 1 - 2 exp(-M m /M^ to 
high redshifts either assuming no evolution in both uq 
and M* (conservative case) or no(z), M*(z) oc (1 + 
z) 145 exp(— z/2.6) (realistic case), which is calibrated 
against Lyman-a observations (refer to for details). 
We adopt a flat ACDM cosmology with Q m = 0.3, 
51a = 0.7, h = 0.7, as = 0.9, the primordial power index 
n = 1, BBKS transfer function[l£| and Peacock-Dodds 
fitting formula for the nonlinear density power spectrum 

The optimal estimator. — Since z > 2 galaxies are 
mainly lensed by matter &t z £ 1 } K = A — B/x(z) 
is an excellent approximation, where A and B are two 
constants and \ is the comoving angular diameter dis- 
tance. Since x( z ) varies slowly at z > 2, one can 
approximate n(x,n) ~ (k) = k((x),?V), where M = 
J2 N[j I J2ij Xj 1 N[j is the effective distance to lens |22J. 
In the limit that Nu ^> 1, Poisson fluctuations become 
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conservative realistic 
n g , {W),{W 2 } n g , {W),{W 2 ) 
123,-0.44,1.1 419,-0.85, 1.5 



76, -0.21,1.2 
40, 0.11,1.6 
13, 0.98,3.7 
3.9, 2.0, 7.7 



290,-0.70, 1.4 
184,-0.49,1.4 
82, 0.1, 1.8 
36, 0.73 2.9 



TABLE I: The predicted number of galaxies n g /l 2 at 2 > 2, 
(W) and (W 2 ) for SKA deep survey with integration time 18 
days/deg 2 . 



stochastic biasing limit [13j where galaxies of different flux 
are uncorrelated with each other. These two cases corre- 
spond to the worst and best cases for the k reconstruc- 
tion, respectively. 

Deterministic biasing. — We first consider the case that 
5 9t ij of different flux bins (but of the same redshift bin) 
are linearly correlated, namely, S g ij = bijSj, where Sj is 
the dark matter density of the j-th redshift bin. As dis- 



Gaussian. The likelihood function of k at an angular cussed above > K can be measured iteratively. Marginal- 
pixel, marginalized over p(<5 9 ,n • ■ • 5 g ,ij), the probability izing over Sj, we obtain 



distribution of 5 g ^j of this angular pixel, is 
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We choose the redshift bin size Az > 0.2 and angular 
pixel size ~ 10 such that 5 3l ij of different redshift bins 
are uncorrelated. For this choice, the matter density dis- 
persion of each redshift bin a m < 0.1. This verifies the 
neglect of high order term S g K in Eq. 1 & 2. 

Since a m < 0.1 and galaxy bias b g is unlikely bigger 
than several [12J, it is reasonable to assume that galax- 
ies are Gaussian distributed. Then p(S g ,u ■ ■ ■) is com- 
pletely determined by the covariance matrix Ci 1 j 1 -i 2 j 2 = 
(^ff,*iii(^)^s,*2j2(^))- SKA can directly and accurately 
measure the correlations of galaxy density fluctuations 
between flux bins, which are the sum of Ci 1 j 1 ^ 2321 cor- 
relations induced by lensing and cross terms. In the in- 
teresting range, Ci 1 j 1 i 2 j 2 dominates. So one can take the 
measured sum as first guess of Ci 1 j 1 -i 2 j 2 . Maximizing L, 
one obtains the optimal estimator k of k. The recon- 
structed k can in turn be applied to subtract the lensing 
contribution in the covariance matrix estimation. This 
can be done iteratively. Since the lensing contribution 
is small, such iteration should be stable and converge 
quickly. 

The properties of high redshift 21cm emitting galax- 
ies are currently poorly known. It is likely that they 
trace the underlying dark matter at some level, and that 
galaxies of different luminosities are correlated to each 
other. We consider this case first, and then the extreme 
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Nij)Nljbij / Nij , Qj = T.i^Wiibii/Nn, u, 
EiWy - N l3 )Nl 3 W l3 /N l3l and T 3 = Ei(^Wy)7% 
N is the mean number of galaxies in each angular pixel. 
(• • • } are weighted by galaxies with the noise from false 
peaks taken into account. 

Maximal stochasticity. — Stochasticity eases the sub- 
traction of the intrinsic clustering signal. In this case, 5 g 
of different bins are uncorrelated. We have 
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FIG. 1: Estimated power of SKA to recover the lensing con- 
vergence map. Left panels and right panels correspond to the 
conservative case and the realistic case. We only use galax- 
ies at z > 2 above a 0.5-cr detection threshold and assume 
unity galaxy bias. In top panels, we fix the integration rate 
as 18 days/deg 2 and plot rms of signal (black lines) and noise 
(shaded regions) as function of smoothing size. Upper and 
lower limits of statistical errors are calculated by Eq. |3]and 
|1J respectively. Bottom panels show the comparison between 
cosmic magnification (fi) and cosmic shear (7) at 10 smooth- 
ing scale. The lower and upper limit of shear measurement 
error are estimated using 5-cr and 10-er cut, respectively and 
adopting mean ellipticity (e 2 ) 1 ^ 2 = 0.54, as expected for disk 
galaxies. SKA is likely capable of reconstructing the k map 
with S/N > 1 at the rate of ~ 4000 deg 2 /year. Cosmic shear 
quickly loses power when decreasing the integration time per 
deg 2 while cosmic magnification is less affected. 



shot noise and the second term is the intrinsic fluctu- 
ation of galaxy number distribution. The conditions 
NijVg,ij ~ O.OlNij -» and N l0 > 1 (for Gaussian- 
ity) can both be satisfied since galaxy bias b g is unlikely 
bigger than several^^- A similar estimator has been de- 
rived by [14}. In two estimators, (W 2 ) and (Wb) arc two 
key ingredients and reflect the key role of flux informa- 
tion. 

Results. — SKA is able to detect n g > 100 arcmin~ 2 
galaxies at z > 2 (table I). For an integration time 
tint = 18 days/deg 2 , a S/N> 2 can be achieved at 
9 ~ 10 (fig. 1). Deep survey configuration detects more 
faint galaxies, which have W — > —2, mimic a constant 
b and thus do not contribute to the signal, due to the 
(W 2 ) - (Wb) 2 /(b 2 ) facotr in Eq. 3. An optimal sur- 
vey configuration should have (Wb) — > 0, which can be 



achieved at t- mt ~ 0.2-1 day/deg 2 (fig. 0. Since n g 
above 0.5-cr decreases much more slowly than t int (for 
example, for the evolution model, decreasing t- lnt from 
180 days/deg 2 to 4 hours/deg 2 , n g only decreases by a 
factor of 9), it is still likely to achieve S/N> 1 at 6 ~ 10 

and scan rate of ~ 4000 deg 2 per year (fig. 2). This will 

1/2 

produce more lensing information (oc S/Nx/ sky ) in a one 
year SKA survey than SNAP(2j| will produce, which will 
cover 1000 deg 2 sky area with S/N~ 2 at smoothing scale 

9 ~ 10'. 

Since the SKA will have ~ 0.3 resolution at z ~ 2, 
it can resolve galaxies and measure cosmic shear. An in- 
trinsic advantage of cosmic magnification measurement 
over cosmic shear measurement is that it does not require 
galaxies to be resolved. Thus, dwarf galaxies which are 
too small and too faint for reliable shear measurement 
still contribute to magnification measurement. Cosmic 
magnification exceeds cosmic shear at integration rate 
< 0.2-10 days/deg 2 (figQ). We note that this compar- 
ison is conservative. We have neglected all systematics 
of shear measurement. For magnification estimation, we 
only select galaxies above a 0.5-cr detection threshold, or 
HI mass above several x 10 8 M Q h~ 2 . There are numer- 
ous galaxies with HI mass ~ 1O 7 M0/i _2 0, which can in 
principle be used to improve the measurement. We do 
not explore its potential in this paper since the luminos- 
ity function at the faint end is unclear. 

Several uncertainties could degrade the signal separa- 
tion. (1) The HI mass function, which is the dominant 
factor, as can be seen from table I and fig. 1. Here 
we further draw the attention on the slope of the HI 
mass function. For an extreme case that a — ► 1 and 
W — > over a large flux range, the signal disappears. 
This effect can be straightforwardly estimated through 
the {W 2 ) and {Wb) terms in Eq. 3 & 4, once the HI 
mass function is measured. Since HI mass function at 
high z is effectively unknown, we postpone the discus- 
sion in this paper. (2) The galaxy bias. For the case 
of deterministic biasing, if b g cx W, flux information is 
no longer useful for the separation and our method effec- 
tively fails. But since b g > 0, as long as the survey is deep 
enough to probe the faint end of galaxies where W < 0, 
b g can not always mimic W and the separation is always 
possible. (3) The galaxy distribution. When b g is big- 
ger than several or smoothing size is smaller than several 
arc-minutes, 5 g is non-Gaussian. In this case, the esti- 
mators described above are no longer optimal. Optimal 
estimators for non-Gaussian galaxy distribution should 
be further investigated. 

Applications. — The reconstructed k map can be ap- 
plied to measure many lensing statistics. For this pur- 
pose, reconstructed k can be noisy because these statis- 
tics generally average over many angular pixels and 
achieve high S/N. Then the optimal estimator derived in 
this paper can be applied to each narrow redshift bins and 
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allows the lensing tomography. (l)The probability den- 
sity function p(n). p(k) as a function of n and smoothing 
angular size can provide independent constraints on 
cosmology. Recently 01 showed that the Wiener filter 
reconstruction of p(n) from noisy convergence map can 
go deep into regions where \k/Ak\ -C 1. We thus ex- 
pect that p{n) can be recovered accurately from SKA. 
(2) Lensing power spectrum and bispectrum. The recon- 
structed k map barely has S/N> 5, so it is consistent to 
neglect < 10% higher order terms: 0{k 2 ) terms and 5 g n 
term neglected in Eq. 1 and k(x) ~ («)• But these terms 
should be taken into account for precision measurement 
of lensing power spectrum and bispectrum ,since their 
statistical errors can reach ~ 1% accuracy |7|. For the 
linear estimator we derived, contributions of these terms 
to the power spectrum and bispectrum can be straight- 
forwardly and robustly predicted. So, there is no need to 
derive a more complicated nonlinear estimator. (3) Clus- 
ter finding and cluster density profile. This is a promising 
approach to break the cluster mass sheet degeneracy. In 
the reconstructed maps, massive clusters at z ~ 0.2 show 
as high peaks with strength k ~ 0.1 and size ~ 10 and 
can be easily identified. These clusters are excellent ob- 
jects to measure the geometry of the universe b y th e tech- 
nique of lensing cross-correlation tomography 16]. Since 
S/N is so high, one can choose smoothing size ~ 1 and 
measure the projected cluster density profile. Exerting 
a prior on cluster density profile, the reconstruction can 
be further improved 17). When k — ► 1, the weak lensing 
condition breaks and Eq. longer holds. By utilizing 
the exact magnification equation, one can develop new 
estimator, in analogy to the reduced shear reconstruc- 
tion 0. We leave this topic for further study. 

We summarize our results. Cosmic magnification is 
statistically more sensitive than cosmic shear because it 
is possible to use the large number of galaxies detected 
at low statistical significance. Intrinsic clustering can be 
subtracted because (1) magnification depends strongly 
on the shape of the luminosity function, which varies sig- 
nificantly, while intrinsic clustering depends weakly on 
the intrinsic luminosity itself and (2) they have differ- 
ent redshift dependence. Cosmic magnification shows 
promise as a complementary technique to map the sta- 
tistically precise distribution of matter, which is not sub- 
ject to most of the systematics of cosmic shear. We have 
worked through the specific numbers for the SKA, but 
the general formalism would also apply to optical spec- 
troscopic or photometric redshift surveys. 
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